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Abstract 

A scattering zipper is a system obtained by concatenation of scattering events with equal 
even number of incoming and out going channels. The associated scattering zipper operator 
is the unitary equivalent of Jacobi matrices with matrix entries and generalizes Blatter- 
Browne and Chalker-Coddington models and CMV matrices. Weyl discs are analyzed and 
used to prove a bijection between the set of semi-infinite scattering zipper operators and 
matrix valued probability measures on the unit circle. Sturm-Liouville oscillation theory is 
developed as a tool to calculate the spectra of finite and periodic scattering zipper operators. 
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1 Scattering zippers 

A scattering zipper describes consecutive scattering events with a fixed number 2L of incoming and 
out-going channels each. It is specified by a sequence {Sn)n=2....,N of unitary scattering matrices 
Sn in the unitary group U(2L) as well as two unitaries U,V & U(L) modeling the boundary 
scattering. The size of the system is supposed to be either even or infinite. Then the scattering 
zipper operator acting on ■ ■ ■ , N}, C^) is defined as 



where the two unitaries Yn and Wat are given by 
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The main hypothesis on each of the scattering matrices Sn is that its upper right entry of size 
L X L is invertible. In the notation of Section [2] below, Sn is in a subset U(2L)i„^ of the unitary 
group. This condition assures that the scattering is effective so that Utv does not decouple into 
a direct sum of two or more parts. The terminology scattering zipper is best understood by 
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Figure 1: Illustration of the first scattering events as well as the boundary scattering. The wave 
functions 0^ = (0^)n>i o-nd ip^ = {'ipn)n>i satisfy 110^ = z(j)^ , which is equivalent to 0^ = z~'^Yil)^ 
and = W0^. The picture corresponds to Proposition [3] below. The boundary conditions 0q and 
ipQ as well as Si can be introduced for convenience, but one can also just keep U , c.f. Section |3| 



looking at Figure 1 illustrating the model. It shows the first scattering events of a semi-infinite 
model N = 00 for which we also drop the indices on U, V and W. It is also possible to consider 
periodic scattering zippers, see Section |8] below for finite operators and Section [9] for infinite ones. 
Furthermore, by placing either of the boundary conditions U and V into Vat, one can consider 
the case of odd A^, but we refrain from doing so. 

The main message of this paper is that scattering zippers are the unitary analogs of Jacobi 
matrices with matrix entries. Here are the structural results supporting this claim: 

• The matrix is not tridiagonal, but five-diagonal. Nevertheless, solutions of the associated 
eigenvalue equation can be calculated by transfer matrices having the same symmetries as 
in the Jacobi case, but no further restrictions (see Section [s]). 

• The matrix element of the resolvent of Uat corresponding to the left boundary site 1 has 
a simple expression in terms of the entries of the transfer matrix, namely it is given by a 
Mobius transformation of the other boundary condition (Theorem [T]). 

• These matrix elements of the resolvent lie on a Weyl surface which is a matrix ball (Theo- 
rems [2] and [3| . 

• All semi-infinite scattering zipper operators U with fixed boundary condition U are in the 
limit point case. Fixing an appropriate gauge for each Sn, the semi-infinite scattering 
zipper operators are in bijection with their spectral measures which are all matrix-valued 
probability measures on the unit circle (Theorem |4]). 

• The eigenvalues of the finite scattering zipper operators can be calculated using matrix 
Priifer phases by Sturm-Liouville type oscillation theory (Theorems [s] and [6]) . This is also 
an efficient tool to calculate the spectrum of infinite periodic scattering zippers (Theorem [T]). 
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For Jacobi matrices with matrix entries, it is well-known how to use transfer matrices, and 
how to calculate the resolvents {e.g. |SB1] . but this is classical). Weyl surfaces are also known 
for Jacobi matrices (see |SB2) which also contains earlier references). The reason why there is 
a simple bijection in Theorem |4] is that all semi-infinite scattering zipper operators are in the 
Weyl limit point case, namely the Weyl surfaces shrink to one point in the large limit. This 
is in strong contrast with Jacobi matrices where there are limit discs which lead to all the issues 
related to the moment problem (see e.g. the beautiful paper by Simon |Siml] ). nevertheless there 
are close connections between probability measures on the real line and Jacobi matrices. Finally, 
oscillation theory of Jacobi matrices with matrix entries was developed in |SBH ISB3| . but again 
the scalar case is well-known. 

Now let us present our personal and without doubt very restricted view on connections of this 
work to the literature. First of all, the scattering zipper itself is a generalization of three well- 
known models, the Blatter-Browne model [BBJ and the Chalker-Coddington network model ^CCj 
of the solid state physics community as well as the CMV matrices |CMV) of the mathematical 
literature. The Blatter-Browne model is scalar, namely L = 1. In the Chalker-Coddington model, 
the aim is to model a higher dimensional lattice of scattering events. In terms of the scattering 
zipper this means that there is supplementary structure in each of the scattering matrices Sn 
and that they are infinite-dimensional and very sparsely filled. Motivated by applications to 
the quantum Hall transitions, the main focus in the Chalker-Coddington model has been on 
random scattering events. First rigorous works on the analysis of the Blatter-Brown model and 
the Chalker-Coddington model on a strip have appeared |BHJt IABJ| . On the other hand, the 
CMV matrices in its matricial version [DPS] only consider scattering blocks of the type 

S = ( ^ (l-ttnO^\ n\ 

where the a„ verify a*Q;„ < 1 which is equivalent to ||q;„|| < 1 (here and below we always use the 
operator norm). These are called the Verblunsky coefficients [DPSJ - There is a huge literature 
on CMV matrices (see |Sim2l IDPS) for a long list of references). They form a subclass of the 
scattering zipper models considered here. Of course, the transfer matrix techniques also apply 
and have very efficiently been used in most works on the subject. However, the Weyl discs at 
finite as presented below seem to have been studied only in the scalar case [GN]. Nevertheless, 
it was possible to prove in [CGZ] IDPS] that every sequence of Verblunsky coefficients corresponds 
to a unique matrix- valued measure on the unit circle (Verblunsky's theorem). Theorem |4] extends 
this theorem in that it exhibits a bijection between all probability measures on the unit circle 
and the semi-infinite scattering zippers with fixed boundary condition U. Finally, it seems that 
oscillation theory for CMV matrices was only developed in the scalar case L = 1 |Sim2[ Theorem 
8.3.3.]. 

Acknowlegements: H. S.-B. wants to thank Mihai Stoiciu for introducing him to the world of 
CMV matrices and for a number of discussions about oscillation theory at a very early stage of 
this work. He also let us know about the reference |BB| . We also acknowlege financial support 
of the DFG. 
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2 Preliminaries on scattering matrices and transfer matrices 

As explained in the introduction, the following subset of the even-dimensional unitary group will 
play a role: 



(3 is an invertible L x L matrix > . (2) 



Proposition 1 /n (|2]), equivalent to the condition that (3 is invertible is either the invertibility 
of 'J or the condition a* a < 1 or the condition 6*6 < 1. Furthermore, one has the representation 

U(2L)i„. = {S{a,U,V) G U(2L) \ a*a<l and U,V e U(L)} , (3) 

where 

S{a,U,V) = ( ^ " , (l-aa*)if/\ _ 
^ ^ \Vil-a*a)-2 -Va*U J ^ ' 

Proof. The equations S*S = 1 = SS* give 

a*a + 7*7 = l, 6*6 + P*/3 = l, 5*7 + = , (5) 
aa* + /3/3* = 1, 66* + 77* = 1 , 7a* + 6(3* = . (6) 

From three of these identities the first claims can be deduced immediately. For ([s]), let us first of 
all note that = 1 — aa* implies that /3* has a unique polar decomposition /3* = U*{1 — aa*)^ 
with some unitary U. Similarly, 7*7 = 1 — a*a shows that 7 = V{1 — aa*)^ for some unitary V. 
But then 6 = -ja* {(3*)-'^ = -V{1 - - aa*)-'2U = ~Va*U. □ 

Recall that the Lorentz group U(L, L) of signature (L, L) is defined to be the set of 2L x 2L 
matrices conserving the form 

^ - [l -°0 ■ (^) 

The following well-known result on the passage from scattering matrices to transfer matrices is 
illustrated in Figure 2. 



" f-f-6(3-^a 6(3-^ 



Proposition 2 The formula 

defines a bijection from U(2L)inv onto \J{L,L). For any vectors 0, 0', ^z^, G , one has the 
equivalence 



Moreover, ifC^' G /orm an inhomogeneity, then 
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Figure 2: Illustration of the scattering event described in Proposition [2| Figure 2(i) is the usual 
graphical representation of a scattering event with incoming amplitudes ip and ip' and outgoing 
amplitudes and 0'. Figures 2(ii) and 2(iii) corresponding to the left and right hand side of 
equation ^ show different representations of the same event as they are used as building blocks 
in Figure 1 and 3 respectively. 



Figure 2 (iii) 



Proof. First of all, note that (p is well-defined because /3 is invertible. Next one readily checks 
that (f{S)*Cip{S) = Chj using the relations ([S]) and ([6]). Moreover, the inverse of Lp is given by 



-1 



A B\\ ( -D-^C D-^ 



^ Dj) \A-BD-^C BD-^ 

Next let us check (|9|, which generalizes ([S]). The upper equation on the l.h.s. is a?/' + (3il)' = + ^ 
which can be rewritten as 

-/3-^a^ + /3-V = i^' -r^i- 

This is already the lower equation of the r.h.s.. Solving it for -0' and replacing in the lower 
equation of the l.h.s. gives the upper equation on the r.h.s.. □ 

Let us conclude this section with a few comments. First of all, one readily checks 



V 0W(l-a*a)-5 (l-a*a)-5a* 
U*) Va(l-a*a)-i (1 - 



Comparing with Q, one way to characterize CMV matrices is therefore to say that its individual 
scattering events give rise to self-adjoint transfer matrices. Of course, this does not imply that 
products of such transfer matrices have the same property. This means that CMV matrices are not 
specified by some symmetry given by a combination of time-reversal, particle-hole or sublattice 
symmetry, and it is not clear to us whether there is a deeper reason to consider scattering events 
of the type |CMV] . On the other hand, it is possible to implement all the above symmetries also in 
scattering zippers, just as for Jacobi matrices with matrix entries (see |SB1| where only even and 
odd time-reversal symmetry is dealt with). Finally, let us briefiy discuss degenerate scattering 
events which don't mix all incoming and outgoing amplitudes. There are many possibilities to 
do this, but only two are relevant for the boundary conditions in the next section. In one only 
depends on ip (left and right in Figure l(i) decoupled), in the other only on ip' (top and bottom 
in Figure l(i) decoupled). Let us focus on the latter. Then there are two unitaries U and V such 
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Figure 3: The scattering zipper of Figure 1 after having transformed each scattering event by (p 
into a transfer matrix multiplication as done in Proposition [3| 



that = Utp' and 



Vtp. In this case, S and ^fiS) are given by 

s = : 1 , ^{s) = 
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3 Solutions and transfer matrices 

In this section, the formal solutions (p^ = {4>n)n>i for the eigenvalue equation lJ(f)^ = z(j)^ at z G C 
will be constructed. Here all the 0^ are L x L matrices and the index n runs to infinity and it is 
formal in the sense that (p^ is typically not square integrable. The construction of (p^ is done such 
that the left boundary condition (at site 1) is satisfied. For finite A^, the solution 0^ in general 
does not satisfy the right boundary condition, namely {1]n4>^)n 7^ ^(plf. 



Proposition 3 Let (pl 

(i) (U0-)„ = z(Pl and 

(ii) (W0^)„ = ^/-^ and {Ytp' 

(iii) For any k > 1, 
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Then the following assertions are equivalent: 
^ = for n>l. 
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where the transfer matrices and initial condition are 



Remark The application (p is used here even if z is not on the unit circle so that z~^S2i is not 
unitary. Of course, 7^^ is in the group U(L, L) only if z is on the unit circle. Also let us point 
out that T^X+i actually independent of z. Moreover, one could also use the initial condition 



and add in (iii) one more equation, namely 
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This results from the discussion at the end of Section |2l o 



Proof of Proposition |3| The equivalence of (i) and (ii) results immediately from U = VW. 
The equivalence of (ii) and (iii) can be checked by successive application of the identity ^ in 
Proposition [2] More precisely, the two equations of (ii) mean that for any k > 1 



V r2k J \<Pk J ' ''^"^ \<p2k+j \r2k+i 

Applying the transformation from Proposition |2] to both of these equations than shows the 
equivalence with (iii). □ 

As usual, the transfer matrices can be nicely iterated. Let us set for n > 

r^{n,k) = r: 

and in order to nicely write the solutions, let us also introduce the notations 
Then the solutions are simply given by 

K = VK-i = r^(n,o)$^, $g = Q . (11) 

4 Resolvents 

As a preparation for the Weyl theory, this section calculates the entries of the resolvents of UAr(y) 
corresponding to the site 1 in terms of the entries of the transfer matrix from to N: 

The first object of study is the Green matrix 

Gl,{V) = 7iI{1]n{V)-z)-'7i, , 

where 7r„ : — )■ ^^({1, . . . , A^}, C^) is the partial isometry onto the nth site and z is in the open 
unit disc D. For self-adjoint matrices, the Green matrix has a Herglotz property, which it is laking 
in the present situation. Actually, for the unitary operator lJ]\f{V) it is more natural to consider 

F^{V) = t7Tl{\JN{V)-z)-\iJN{V) + z)7r,. 

One readily checks that this analytic function 2; G D F^{V) G Mat(L,C) has a positive imag- 
inary part t{F{z)* — F{z)) > 0. It also satisfies -F(O) = il and is thus a so-called Caratheodory 
function (up to the factor i). If one also takes the Cayley transform of the domain D to the upper 
half-plane H, one has again a Herglotz function. A little more information on this is resembled 
in Appendix B. 
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Theorem 1 The matrix Cf^ — VAf^ is invertible and 

lies in the Siegel disc = {Z G Mat(L,C) | Z*Z < 1}. Moreover, the resolvent matrix F^{V) 
is given by 

Fn{V) = \ (EUV) + 1) (EUV) - 1)-' , (12) 

and the Green matrix by 

Gl^iV) = EUV) (1 - EUV))-' . (13) 

Let us note that the formulas for Elf{V) and F^{V) can also be written using the (inverse) 
matrix Mobius transformation (see Appendix A): 

EUV) = V* : r'{N,0) , F^{V) = C* ■ E%{y) , (14) 

where the Cayley transformation is the 2L x 2L unitary matrix defined by 

C = ^(^ -'^ 

72 VI ^1 



Based on (14) and the results of Appendix A, one obtains further identities, e.g. 

E%{V) = r'iN,0)-' -V* , F^{V) = {C*r'iN,0)-^) -V* . (15) 

For the proof of Theorem [1} one needs a number of lemmata which will also be useful later on for 
other purposes. 

Lemma 1 Let = Lp(z~^S) for some S E Uin^(2L) and 2 G D, and set 

^ B\ _ f{\z\-'-l)A*A {{z)-'-z)A*B 



C zDj ' \{z'^-z)B*A {l-\z\^){B*B + 1] 

Then > ^— ^ and 

iry = c + v . (16) 

Moreover, for any Z G Dl, the Mobius transformation (T^)~^ ■ Z is well-defined and lies in D/,. 
Furthermore, only for 2; G D, 

(T')-^ -Z e ^L- 

Proof. With the notations of the lemma, the matrix T = ^{S) satisfies the defining equation 
T*CT = £ of the group U(L, L) so that 

^^^•£(;^^U£. (17) 



CD \C D 
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Thus a direct computation using the identities contained in (17) leads to (16). Note that alter- 
natively can be expressed in terms of C and D only. It remains to show the positivity of V^. 
Let (J) G C^^ be with, say, ||^|| > Then 



= {\z\-^ - 1) (f)*A*A(j) + {{z)-^ - z) <j)*A*Bi) 

+ [z-^ - z) i)*B*A(j) + (1 - \z\^) i>*B*Bi> + (1 - \z\^) i>*i) 
> {\z\~^ - 1) (f)*A*A(j) - 2 |(^)-^ - z\ {ip* B* B^)'^ 



+ (1 - \z\^) ^Ij*B*B^Ij + (1 - \z\^) ^*ip 



> 



{\z\-^ - ly {4>*A*A4>y^ - (1 - \z\^Y {i)*B*Bi)Y^ + (1 - \z\^) i/j*^ 



1p J \1p 



> 



where in the last inequality the bound 
same manner. 

As to the last claims, let us begin by noting 



(r 



was used. The case 



zA* -C* 
-B* ^D* 



< 



is dealt in the 



(18) 



Thus we need to check the invertibility of —B*Z + j;D* for Z G and z e3. For that purpose 
let us use again (17) so that, in particular, D*D — B*B = 1. Thus D*D > and D is invertible. 
Also {BD-y-BD-^ = 1 - {D-y-D-^ < 1 from which follows \\BD-^\\ < 1. Therefore 

-B*Z + - D* = - D* (l-z (BD-^yz) , 

z 



z 



is indeed invertible. Thus (T^)"^ ■ Z = {zA*Z - C*){-B*Z + z-^D*y^ is well-defined and 



But using (16) 



l] ((rr^)*/:(rrMi 



z 



1 



z 



{-B*Z + z-^ D*) . 



(19) 



< 



z^ 



Z*Z - 1 < 



which together with (19) shows the last claim. 



□ 



Corollary 1 With the positive matrix V2}. defined in terms of entries of T^k in Lemma [TJ 



N/2 



T'{N,oycr{N,o) = c + Y.O'ir---{%i_,yv^,,%i_,---r,^. 
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Proof. This follows by iterating Lemma [T] and using (7^fc_i)*'C7^|._^ = £. □ 



Lemma 2 Let a,b be L x L matrices and set $ = (^) . If ^*C^ < 0, then b is invertible. On the 
other hand, if > 0, then a is invertible. 

Proof. < implies that a*a — b*b < 0. If there is a vector v such that bv = 0, then 

v*a*av < which is impossible and therefore b is invertible. If a*a — b*b > 0, the same argument 
implies invertibility of a. □ 

Lemma 3 Let ^ and $ be two 2L x L matrices of maximal rank satisfying \I'*£\[' = and either 
> or < 0. Then '^*C^ is invertible. 

Proof. The claimed invertibility does not depend on normalization so that we may assume that 
\E'*\I' = 1 and $*$ = 1. Now the fact that is £-Lagrangian implies 

+ = 1 . 

Applying \&* and $ to the left and right of this equation shows 

Thus let us argue that <|)*\E'\E'*$ < 1 because this then shows that $*£\E'\E'*£$ and thus also 
\E'*£$ is invertible. Now <|)*\i>\['*<|) is given by the squares of the scalar products of the vectors 
in the planes spanned by $ and \I' and its eigenvalues are thus the squares of the cosines of the 
principal angles between these planes. There is an eigenvalue 1 if and only if one angle vanishes 
and therefore if and only if $ and have a direction in common. This would mean that there 
are non- vanishing vectors v,w & such that "^v = $w, which is incompatible with \E'*£\E' = 
and > (or < 0). □ 

Lemma 4 Suppose \z\ < 1. Given ^ = {C,k)k=i,...,N with C,k ^ Mat(L,C), the solution of the 
equation (Uat — z)(j) = ^ is, for even n, given by 

where o-nd 62k o-f^^ the entries of S2k o-nd (pi has to be chosen in the unique manner such that 



ipN given by (|20|) satisfies 

/1\ ^'^ f- - 

-^k 




Proof. As above, the equation (Uat — z)(j) = is solved using an auxiliary vector ip satisfying 
Yip = z(p + ^ and W0 = ip. This is equivalent that, for k = 1, . . . , y and respectively k = 



1 ^ - 1 

J., . . . , X 



' 2 

_ T 

' 2 
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together with the condition ipi = U(f)i and the constraint that ipN = V(f)N stemming from the 
other boundary condition. Each of the two equations in ( [22| ) (the first one divided by z) is 
transformed using Proposition |2| 

^2 J " '^'\ct^2k-l)^\ -P^k )\i2k) ' \ct^2k+l) ~ ''''''\^2k 

Iterating and replacing ipi = Ucpi gives, for even n, 



while for odd n the entries on the l.h.s. are simply exchanged (c/. the definition (10) of $^). 
Now one, moreover, has to satisfy the constraint V(j)N = 'ipN- For that purpose, one takes the last 
equation for n = and extracts ^/'jy which is then set equal to V(f)N. This leads to equation (21 ) 
which can indeed be solved for 0i because the matrix 

o) 

is invertible by Lemma |2j Indeed, Corollary [l] shows that $ = TJ^---7^^(^) satisfies > 

which is the hypothesis in Lemma |2j □ 

Proof of Theorem [l] By the last claim of Lemma [l) (7^_i)-^ ■ ((7^)"^ ■ 1^*) = (7^7^-i)"^ ■ V* 
is well-defined and lies in the Siegel disc D^. Iterating this shows that 

{%'r,T' ■ (■ ■ ■ ((r^r^_i)-' ■¥*)■■■) = r'{N, oy' -v* = v*-. r'{N, o) , 

exists and lies in Dj;,. This shows the first claim. Let us note that the invertibility V*C^ — = 
V*{Cf^ — VA^j^) also follows from the identity 

because $ = 7^ ■ ■ ■ 7^^ (J) satisfies > as shows Corollary [l| so that Lemma [s] applies. 

The Green matrix Glf{V) is the component (pi of the solution of {IJn{V) — z)(j) = C, with 
inhomogeneity = ^i.fcl- Hence by Lemma |4| 



l^ n-z n-z f A^ ^z nr\ , fn-zn-z\-l ( — ^ ^<^2/32 ^ ^ A 



and 



Now as (7^^) = C{T^ one checks that using the identities in U(2L)i„ 

{%'riT 
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Therefore the above two equations become 

<Pn = (A^ + Blf) G'^iV) + z-^Bl, , V^N = {C^N + D%) GUV) + z-^D'^ . 
Now these two equations have to be solved for Glf{V). For that purpose, the invertibihty of 

is needed. It foHows from Lemma [s] because ) is £-Lagrangian and ^ = T]^ ■ ■ ■ l~i satisfies 
> as already argued in the proof of Lemma |4j Therefore 

G%iV) = ^ [iD% - VBl,) + {C% - VA^Y'iVB^ - D^,) . 

Now the formula for G^j^{V) follows. 

Finally let us calculate F^{V). We start from 

iF^V) = {tiUVj^- z)-^n2n;S2n^ + n^il]J^ - z)-\^ S2ni)U + z 71,(1] n - z)-^ 

= [ttHVm - zy^TT^ 72 + Gr,{V) a2)U + z G^V). (24) 

Let us denote Gi 2 = ttKVn — z)^^tt2, thus we can calculate Gi^2 using the same procedure than 
with Gj\f{V). Notably, G12 is the equal to the component 0i of the matrix- valued solution of 
(Uat — z)(j) = with inhomogeneity = 52,fcl- Hence 



and 



Using again (T;^)-! = CCJ^ ')*C and the identities in U(2L)i„„ one checks that 



Therefore the above two equations become 
and 

V<Pn = {C'n + D%)G,,2 - C'nU*12' - z-^D'^a2i2' ■ 

Solving yields 

Gi,2 = [{D^^-VB^U + {C^^-VA%)Y\{D^^-VB^^){-z-'a2i2')-{C^N-VAUU*^2'] 



Replacing this and G^iy) in the form (13) into (24) leads to the desired formula for F^iy). □ 
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5 Weyl theory 



For z G D, let us set 

We consider Qj^ as quadratic form on and are particularly interested in its maximally isotropic 
subspaces, also called Lagrangian subspaces. It will shortly be shown that the signature of is 
(L, L) so that the dimension of these Lagrangian subspaces is L. By definition the Weyl surface 
is then the image of (the inverse Cayley transform of) these subspaces under an adequate chart 
on the Grassmannian called the stereographic projection. Recall that the Grassmannian of all 
L-dimensional subspaces of C^^ can be introduced as set of equivalent classes of 2L x L matrices of 
maximal rank with respect to the equivalence relation $ ~ v]/ <^=r^ $ = v|/c for some c G G1(L, C). 
Elements of will be denoted by [^]^. Let us consider the subset Q'l" C G/, of subspaces 
represented by some $ = with an invertible L x L matrix b. This set is the domain of the 
stereographic projection tt : G^" — ?■ Mat(L,C) defined by 

nii^U = ah-\ $ = f ?V (25) 



Definition 1 The Weyl surface is defined by 

dW'^ = 7r({[C*<l>]^ G I $ isotropic for Q%}) , 

and the closed Weyl disc by 

W% = n{{[C*^U e Ol I < 0}) . 

Of course, we have to check below that this is well-defined, namely that [C*$]^ is in the 
domain of vr for all Lagrangian subspaces of Q^. Then the next aim will be to show that the 
Weyl surface is the surface of a matrix ball and this will ultimately allow to prove estimates on 
the dependence of Glf{V) and F^{V) on the boundary condition V G U(L). 

Proposition 4 The quadratic form has the following properties: 

(i) Qk+i = Qk ^nd Ql, = Ql,_, + {r,r ■ ■ ■ {U-iyVk%U ■■■r,\ Furthermore 

N/2 

= £ + Y^iT.r ■ ■ ■ {%ur vk %u ■■■v 

k=l 

where is positive matrix given in term of entries of 7^^ defined in Lemma [T]. 

(ii) {Q%y' = cQi-'c 

(iii) Qk > Qlk^i 

(iv) signature((5|^) = {L, L) 

Proof, (i) is just a reformulation of Lemma [l] and Corollary [l| (ii) results from (18). (iii) is a 
consequence of (i) and > 0. Finally, from the definition of Q^, one deduces that 

signature(Q^) = signature(£) = (L, L) , 

which is (iv). □ 
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Proposition 5 Let Ef^iV) be given as in Theorem [Tj // $ is Lagrangian for Q}^, then th 
a unique unitary V and an invertible matrix b with 



ere is 



Moreover, [C*$]^ is in the domain ofn. 



Proof. Let $ = (^) satisfy = 0. This implies that < by Proposition [ij^i) and 

thus by Lemma [2 that b is invertible. As (^) is isotropic if and only if (^''^ )6 is isotropic, we may 
assume that 6=1. Now $*Q|^r$ = is equivalent to 



Using polar decomposition, it thus follows that there exists a unique unitary V such that 



Dn - VB^^ 



As VA^f — Cf^ is invertible by Theorem [T| it follows indeed that a = Ef^{V). Furthermore, 



EUV) + 1 
z{El,{V) - 1; 



(26) 



By Theorem [T| ^^(V) is in the Siegel disc so that El^{V) — 1 is invertible and therefore [C*$]^ 
is in the domain of vr. □ 

Theorem 2 One has 

= {F^{V)\VeV{L)} = {t{l + 2zG%{V)) \ V eV{L)} . 



Proof. The first equality follows immediately from (26) and (12) in Theorem [Tj The second one 
then follows by combining (12) and (13). □ 



The next aim is to analyze the geometry of the Weyl surface. This can be done in complete 
analogy with |SB2] if one works with the Cayley transform of the quadratic form: 



N/2 



Q-^ = C*Q%C = tJ + J2 C*T'{2k-l,oyv^i,T%2k-l,0)C 



(27) 



fc=i 



where 



J 



C*CC 



-1 

1 



Definition 2 The radial and central operators are defined by 

-1 



Qn 



S 



N 



R 



N 



Q 
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Proposition 6 Let z G D. Then Rf^ is well-defined, positive and decreasing in N. Also Rf^ 
well-defined, but negative. Moreover: 

( QZ \* QZ~^ 

Wn) — '-'N ) 



IS 



(S'^)*(i?^) + R 



N 



Proof. All claims on Rf^ follow by taking the matrix element of (27) because (J) J7'(J) =0 
and the map r G (0, oo) h- )■ —r~^ is operator monotone. For the second claim, one has to adapt 
Lemma [T] and Corollary [T] to the case \z\ > 1, but then the proof is identical. Next let us note 
that the Cayley transform of Proposition |4](ii) reads 

J = Q'n'JQn ■ 

Using :r =(;)Q*- (J) g)*, this yields 



J 



N 



Qn - Q 



N 



Q 



N ■ 



1/ \o, 

The upper left and upper right entries of this equation give 

= -{Rlr )~^S^ (Rn)'^ + (R-N )~^('S'^)*(-RAr)~^ 
-1 = {R% r^S% {R^y^S^^ - {R% y^i-A ^atL 



These two equations lead to the remaining two claims. 



□ 



Theorem 3 The Weyl discs are strictly nested matrix discs, namely one has QU^ C 2n^_]^ and 
(92n^_i n Wj^+i = 0, as well as 



5^+ {R%r^W{-R 



■N 



W G U(L) 



Proof. The first claim is an immediate corollary of Proposition 4]^iii). Let = (f) be a 
Q^— Lagrangian plane. This is equivalent to 



= iE+ lYiR^r'iE +1) + ^(E- ir{s^^r{R^^r\E + 1) 

-z{E + inR%)-'SUE-l) + {E-ir{{S%r{R%)-'S^^ + Rl-'){E-1), 
where in the second equality Definition |2] and Proposition |6] were used. Rewriting gives 

(^{E + l)-^SUE-^))\R%)-'[{E + l)+^SUE-l)) = - {E - 1)* R'^-'iE - 1) 
By Proposition [5] and Theorem [TJ G D/^ so that — 1 is invertible. Hence 



tiE + l)iE - 1)-' - St,) {W^)-H-i{E + 1){E - - S 



'AT 



-R 



N 
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Therefore there exists a unique unitary W such that 



I 

As the l.h.s. is precisely 7r([C*$£;]^), this concludes the proof. □ 
The following is now an immediate consequence. 

Corollary 2 For any boundary condition V and 2; G D, there is a unitary W G U(L) such that 

fn{v) = s't, + {R^^y^wi-Ri-y^ , 

and 

G%{V) = l-(l-^S%- ^{R^^y^W{-Rl"y^) . 

By now the beautiful theory of Weyl discs is complete. Let us finally come to its main 
analytical application, namely the control of the dependence of the Green matrix Gf^{V) and the 
resolvent F^{V) (both matrix elements of the resolvent at the boundary site 1) on the boundary 
condition V on the (other) boundary site A^. Let us start by noting the following. 

Corollary 3 For any V, V G U(L), 

This becomes useful in combination with the following result. 
Proposition 7 One has 

~i 18 



^" ' " ^ "/ - (1 - \z\^y N 
Proof. Using (27) and the bound from Lemma [l| 



(0)^^(0) = l^i^^ ^^^'y--<^2i-irv^,T,u---v 



4 

k=l 



Now ($^)*$f = 2 1 and for > 2 



mm > mr^n = Q) q^Q) > Q) ivrv^r.^l^l] > (i-k^i. 

(This is rough, but sufficient for our purposes.) Therefore 



from which the upper bound on R^ follows. Using a version of Lemma [l] for 2; ^ D, the bound 
on i?^ is shown similarly. □ 
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6 Spectral measures of semi- infinite scattering zippers 



Let now U be a semi-infinite scattering zipper operator associated to a sequence {Sn)n>2 of 
scattering matrices Sn € U{2L)i^^ and a boundary condition U G U(L). Then it follows from 
Proposition [7] combined with Corollary |3] that the limit 

= hm F^iV) , 

exists and is independent of the choice of V. In the terminology of Weyl theory, a semi-infinite 
scattering zipper operator is always in the limit point case. Moreover, the convergence is uniform 
in z on compact subsets of D. Therefore, 2; G D i— )■ is analytic. As '^m{F^{V)) > and 
F^[V) = tl it follows that also 53m(F^) > and F^ = il. Therefore the Riesz-Herglotz 
representation theorem recalled in Appendix B can be applied to show that there is a unique 
matrix-valued probability measure /i on §^ such that 

F^ = I ! f^ido 1^ . 

This measure is called the spectral measure of U. It dominates all other spectral measure (because 
the range of tti is a cyclic subspace for U). Resuming, there is a map that associates to every 
semi-infinite scattering zipper operator a matrix- valued probability measure n = '(9(U) on 
There is still some gauge freedom though allowing several operators U to have the same spectral 
measure, as already pointed out in [BB^ iBHJj . This is dealt with in the following result. 

Theorem 4 Let Ai{U) denote the set of all semi-infinite scattering zipper operators with left 
boundary condition U G U(L) and scattering matrices {Sn)n>2 given by = S{an, Un, Vn) where 
< 1 and Un,Vn G SU(L). Then the map i} establishes a bisection between M.{U) and the 
matrix-valued probability measures on 

For the proof, it remains to construct an inverse to namely to construct a scattering zipper 
with boundary condition U from a given matrix- valued probability measure yU on S^, such that /i is 
again its spectral measure. Adapting the approach in [ DPSj . this will be done by a Gram-Schmidt 
procedure associated an adequate scalar product ( • , . ) (with values in the Lx L matrices) and an 
adequate basis for the functions on §^ (this will be given by Laurent polynomials in an adequate 
order as suggested by |CMV| ). Let us define 

{f,g) = [ g{z) f,{dz) f*{z) , (28) 

where / and g are matrix- valued functions. At first sight, there seems to be something wrong with 
this definition. Indeed, one can define another scalar product by exchanging / and g on the r.h.s.. 
When the Gram-Schmidt procedure below is done w.r.t. this other scalar product, one obtains 
different orthonormal polynomials, and only their so-called Szego transformation |DPS] will lead 
to the polynomials (j)^ and ip^ used in SectionJS Somewhat anticipating, already the notations 
(j)^ and will be chosen already now. Using (128) allows to avoid using the Szego transformation 



and also simplifies the calculations below. The product (28) is left matrix-linear in the second 
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argument and anti-linear in the first, namely it satisfies for matrix valued functions /, g and h 
on as well as a matrix a that 

{f,g + ah) = {f,g)+a{f,h) , {f + ag,h) = {f,h) + {g,h) a* . 

Now two sequences 0^ = {(f)^)n>i and ip^ = {'ipn)n>i of orthonormal families of matrix-valued 
Laurent polynomials in z will be constructed by the Gram-Schmidt algorithm w.r.t. (28). For 
0^, one orthonormalizes the sequence {1, z~^l, z^l, z~^l, z"^!, . . .}. Then (pl = 1 and the nth 
element of the resulting orthonomal sequence is 0^. Simarly, tp^ is obtained by orthonormalizing 
{U, zl, z~^l, z'^1, z~'^l, . . .}. Hence ipl = U . The orthonormality relations read: 



Let us note that the (p^ satisfy 

<Pln = t^2n + V{z) , 4>2n+l = '^2n+l + q{z) , (29) 

where the leading coefficient «;„ are invertible matrices, p is a matrix polynomial in the span of 
2;""^, . . . , z~'^~^^ and q a polynomial in the span of . . . , Similarly, 

^2n = K,2n Z"" + p{z) , tp^^+l = K2n+1 Z~ + q{z) , (30) 

with invertible /t„ and polynomials p and q in the span of ^, . . . , z and z"^,.. . ,z 
respectively. Now let us also define 

Pn — '^n-l('^n) ; Pn — l^n-l{,l^n) 



Lemma 5 The orthonormal polynomials (p^ and ip^ have the following properties. 
(i) There exist L x L matrices Q^n+ii ^"in+ii ^\n+2i ^2n+2' such that 

- P2n+102n+l = "L+l'/'L . '/'L " P2n+1^2n+l = C^ln+1^2n ^ (31) 



and 

^^"VL-I -P2n02n = «L02n-l , #2n-l " P2n ^2n = «2n ^2n-l • (32) 

(ii) One has 

1 ('2 \* 3 / 4\* 

'^2n+l ~ \'^2n+l) 5 '^2n ~ l0^2nJ • 

In the following, we will simply denote «2n+i '^2n+i (^iT'd ot\^ by a2n- 

(iii) There are unique unitaries Un and Vn with unit determinant such that 

pn = (1 - a„a*)5[/„ , Pn = {l-alan)^V*. (33) 

(iv) One has a*a„ < 1. 
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Proof. All four relations in (31) and (32) are dealt in the same manner, so let us focus on 
the first one. By (29) and (30), — p2n+i4>2n+i ^^e span of z""'^, . . . , z~". Moreover, 

this polynomial is orthonormal to z"~^l, . . . as ip2n ^i-nd 02n+i by construction. This 

implies that -i/^g^ — P2n+i4>2n+i ^ ^^^^ multiple of 
For item (ii), let us first calculate as follows: 



a 



1 

2n+l 



an. 



2n+l\r2n) V2n 



02n) = (02n, "L+l^D = (02n, ^2n " P2n+102n+l) = (0L,V^2n) 



On the other hand. 



a 



2n+l 



Thus af^^i 



«L+l(V'2n,V'2n) = (^2n, «2n+1^2n) = (V^2n, 02n " P2n+1^2n+l) = (^2n, 02n) • 
1 



2n+lJ 



The other equality is checked in a similar manner. 



Next let us check the first identity of (33) for odd index: 

= (P2n+10L+1 + a2n+102n' P2n+102n+l + a2n+10L) 
= (P2n+102n+l)P2n+102n+l) + (tt2n+102n) '^2n+102n) 
P2n+1 (02n+l 5 */'2n+l)P2n+l 

+ a2n+l(02n>02n) "2n+l 

= P2n+lP2n+l + tt2n+l"2n+l • 

Thus, there exists a unique unitary of unit determinant denoted U2n+i such that 

P2n+1 = (1 ~ (^2n+lCX*2n+l)'^ f^2n+l • 

The other cases are dealt with in the same way. As each p„ and p„ is invertible, the identities 
(33) also imply (iv). □ 

Proof of Theorem|4j From the density of Laurent polynomials, it follows that (0^)n>i and ('?/'^)n>i 
are orthonormal basis of the square integrable matrix- valued functions on w.r.t. ( • | • )• Thus 
any L x L matrix valued function / can be expanded as follows: 



n=l 



oo 

E 

n=l 



(34) 



Now let us define the matrix entries of semi-infinite matrix-valued operators by 



Then it follows from (34) that 

oo oo oo 



u 



n,m 



k,m 



k=l k=l k=l 

It remains to show that V and W defined above have the same structure as their homonyms from 



Section |3| For n > 1 it follows from (31) that 

a2n+102n + P2n+102n+l , '^2n+l = {P2n+l)* <p2n ~ ^2n+lO*2n+lU 2ri+l<p2n+l 



2n 
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where the identity (p2n+i) ^Oi2n+iP2n+i = V2n+ia2n+if^2n+i was used. Recalhng the notation Q, 
it follows that 

iz^" ) ^ S{a2n+l,U2n+l,V2n+l)i .^'^^ 
K'r2n+1/ \V2n+l/ 

Together with Wi^i = {(j)\^ipl) = U, one obtains 

W = f/e0 5(a2fc+i,f/2fc+i, W) . 

k>l 



Similarly, (32) implies 

2^" VL-I = «2n02n-l + P2n02n , Z'^i'ln = (P2n)>2n-1 " '^2na2nf^2n02n- 

Thus, for n > 1, 



] = Z 6(a2n, U2n, V2nj 



and one obtains 

V = S{a2kU2k,V2k) . 

k>l 

In conclusion, U is a semi-infinite scattering zipper operator. □ 



7 Intersection theory and oscillation theorem 



The transfer matrices allow to calculate formal solutions of Uat^^ = Z(f)^ by equation (11). As 
already pointed out in Section [3| the ^-component of $^ does not yet lead to a matrix-valued 
solution of Utv^^ = z(j)^ because there is a supplementary constraint on the two components of 
Indeed, because is even, the last equation of Wtp^ = (j)^ is V(f)]^ = ip^. For each vector 
V G satisfying Wip^v = (p^v one finds an eigenvector of Uat. Let us reformulate this in terms 
of the dimension of an intersection of two L-dimensional subspaces of C'^^: 

multiplicity of z as eig envalueofU^v = dim($^C^ n ^yC^) , (35) 

where 

Now the intersection of the planes $^C^ and \E'yC^ can be conveniently calculated using the fact 
that both planes are £-Lagrangian for z G namely they both satisfy = 0. In fact, the 

initial condition $q as well as \E^y are £-Lagrangian as shows a direct calculation. As the transfer 
matrices for ^ G all conserve the form C, also $^ is £-Lagrangian for all n. Before analyzing 
the intersection of two £-Lagrangian planes, let us study the set C of Lagrangian planes, 
namely those [$]^ G Gl satisfying = 0. 

Proposition 8 hi G G^" and tt : ^ U(L) is a bijection. Every [$]^ G has a representa- 
tion of the form $ = (^) with U G U(L). 
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Proof. Let $ = (^) be £-La grangian. Then ker(a) fl ker(6) = {0} because otherwise $ would 
not be of rank L. Moreover, a*a = b*b so that both a and b are invertible. Also U = ab~^ G U(L) 
and $ = b. From this all claims follow. □ 

Now follows a general result about the intersection of two £-Lagrangian planes. 

Proposition 9 Let $ and \E' be C-Lagrangian frames and set W = 7r([$]^)*7r([^]^). Then 

dim($C^ n \I^C^) = dim(Ker($*£\E') ) = multiplicity of 1 as eigenvalue of . 

Proof. Let us begin with the inequality < of the first equality. Suppose there are two 2L x p 
matrices v,w of rank p such that $f = ^w. Then $*£\E'w = = so that the kernel 

of $*£\E' is at least of dimension p. Inversely, given a 2L x p matrix w of rank p such that 
$*£\E'w = 0, one deduces that (£$)*\E'w = 0. As the column vectors of $ and are orthogonal 
and span C^, it follows that the column vectors of "^w lie in the span of $, that is, there exists an 
2L X p matrix v of rank p such that \E'w = $f . This shows the other inequality and hence proves 
the first equality of the lemma. For the second, let us first note that the dimension of the kernel 
of $*£\E' does not depend on the choice of the representatives. Choosing the representatives 
$ = (''(I*!-)) and ^ = (''(["J^l-)) then shows the second equality. □ 

For 2 G let us define G U(L) to be the unitary associated by Proposition^ to $ = $|^r 



and \E' = Taking into account the explicit form of \E'v' as well as (10), one finds 



= rNi<p%r'v , (36) 

These unitaries are the analogs of matrix-valued Priifer phases and can conveniently be calculated 
by iterated Mobius transformations with the transfer matrices. Now the main facts of oscillation 
theory can be stated. 

Theorem 5 Let N >2 be even. For z eS^ , one has 

multiplicity of 1 as eigenvalue of = multiphcity of z as eigenvalue of Uat . (37) 
Furthermore, setting z = e*^, 

- (W^ydeW^ > , (38) 
t 

so that all eigenvalues of rotate in the positive sense as a function ofOES^. 



Proof. The equality (37) follows immediately from (35), the definition of and Proposition^ 



Further, let us note that = ■?/'^(0^) = ((V'^)*) ^{4>lf)*V. Thus one calculates 

-^{W^TdeW^ = ^l^*((0^)"^)*[(^^)*9,^^-(0^)*9,0^](0^)-V 
= V*{i<Pl,)-'yi<^%ntC){de^l,) (0^)-V . 



It is therefore sufficient to check the positivity of ($^)*(z£)c?6i$|^. By ([TT]) and {T^)* CT^ = 
one now finds 



N 



.n=2 
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For odd n, the summands vanish because is independent of z and thus also 9. We shall now 
show that for every even n there is a strictly positive definite contribution. For that purpose, let 
us drop the index n on and Sn- Now = ^{zS) is given by 

^ (zA B 
\C zD 

where A, B, C and D are the coefficients of = ^{S) G U(L, L). Then one verifies 

It remains to check that this matrix is positive. On first sight, it is not even hermitian, but 
actually the defining relations of U(L, L) are 

A*A-C*C = 1, D*D-B*B = 1, A*B = C*D . 

The last one shows that the above matrix is indeed hermitian. Now let (^) G C^^ be with, say, 
> U\\. Then 



A* A zC*D 
^P) \zB*A D*D J U 



'A*A(j) + z(t)*A*Bip + z^Ij*B*A^ + iIj*B*B^Ij + ip*ilj 



> 

1 

> 

~ 2 v^y 



where in the second step we used the Cauchy-Schwarz inequality. This completes the proof. □ 

Remark As Uiv has exactly NL eigenvalues, the rotation number of 2; G §^ )• is equal to 
NL. This can also be shown independently by calculating the Maslov index as in |SBlj . o 



8 Oscillation theory for finite periodic scattering zipper 



It is possible to associate to a sequence {Sn)n=i,...,N of scattering matrices Sn G U(2L)inv with 
still even, a periodic scattering zipper operator U^" = V^"W^'' by setting 



(82 



■n rper 
^ N 



\ 



JPCT 



where 



Sn/ 
Si 



71 



S, 



Sn-1 



aij 



"1 /3i 
7i ^1 
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The aim of this section is to calculate the spectrum of This parallels the calculations in 

|SB3] and is useful for the calculation of the spectrum of infinite periodic scattering zippers, 
as explained in Section |9} The solutions 0^ G ^^({1, • • • (g) of the eigenvalue equation 

U^''0^ = z(f)^ can again be constructed with the transfer matrices. Indeed, every eigenvector of 
T^{N,0) to the eigenvalue 1 allows to construct a periodic eigenvector so that 

multiplicity of z as eigenvalue of U^' = multiplicity of 1 as eigenvalue of T^{N, 0) . 

Hence one needs to find those z G for which 1 is eigenvalue of T^(iV, 0). By considering the 
fact that the graph of T^{N, 0) is also Lagrangian w.r.t. an adequate quadratic form, this can be 
done by means of intersection theory in a similar manner as in the previous section. 

Several notations need to be introduced. Let us associate to T the 4L x 4L matrix T = l2L©7~ 
where the © denotes the checker board sum given by 



(39) 



This allows to define the quadratic form C = C® C. For z G §^ the matrices conserve the 
quadratic form C Therefore also £-Lagrangian planes are mapped onto £-Lagrangian planes. 
By Proposition [sj the stereographic projection maps the set L2L of >C-Lagrangian planes is diffeo- 
mophically to U(2L). Let us denote the stereographic projection by 7? : L2L — t- U(2L). 

Proposition 10 To a given T G U(L,L), let us associate the unitary 

W = 7r([f$o]~)* v?([$o]~) G U(2L), 
where the C-Lagrangian plane is given by the AL x 2L matrix 





B \ . 






it 




D 1 c' 


D' 



( ^ 





B 





\ 





A 





B' 




C 





D 







\ 


C 





D' 


I 



(40) 



/ 


1 \ 


1 





1 





^0 


1/ 



(41) 



Then 



geometric multiplicity of 1 as eigenvalue of T = multiplicity of 1 as eigenvalue of W . 

Proof. First of all, it can readily be checked that \1/q£ \1/o = 0. Now let us suppose that this 
frame and the Lagrangian frame T ^^o have a non-trivial intersection. This means that there exist 
vectors v^w^v'^w' G C'^ such that such ^l'o(^) = '^^o(^/)- The first and third line of this vector 
equality imply w = w' and v = v', the other two that T(^) = (^). This shows 



geometric multiplicity of 1 as eig envalueof r = dim(r^oC^^ n ^0 C^'^) 
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But now Proposition |9] can be applied to calculate the r.h.s. and completes the proof. 
It is now natural to introduce the following £-Lagrangian frames: 

$^ = $^ , $^ = $n n > 1 

'n ^ n-1 ' ^0 ^0 5 ^ 



□ 



(42) 



with \E'o as in (41). Associated are then the unitaries 

Theorem 6 The multiplicity of z = e*^ G S"*^ as eigenvalues of U^'' is equal to the multiplicity of 
1 as eigenvalue ofW^. Moreover, 

- [W^ydeW^ > 0, 

so that the eigenvalues of rotate around the unit circle in the positive sense and with non- 
vanishing speed as function of 6. 



Proof. The first claim follows directly from Proposition [TOj For the proof of the second one, let 
us denote the upper and lower entries of \E'^ by ip^ and tpt- These are 2L x 2L matrices such 
that 7r([$^j^) = = ((^l)-i)*(^/'^)*. Then 

Thus it is sufficient to check positive definiteness of 



N 



From the product rule follows that 



N 



N 



n=l \/=n+l 



This implies that 



N 



'n-1 



n=l 



As 



and the matrices do not mix first and third columns and lines with the rest, it follows from 
evaluation in the state \E'o that 



N /n-l 



fn-1 



n=l \l=l 



But positivity of {T^)* (iC) d^T^ was already checked in the proof of Theorem [5j This completes 
the proof. □ 
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9 Spectrum of infinite periodic scattering zippers 

In this section, we consider a two-sided infinite scattering zipper W defined on e^{ZX^) which 
are iV-periodic where again is even. It is specified by a sequence of scattering matrices Sn € 
U(2L)i„^ satisfying Sn = Sn+N for all n G Z. One can partially diagonalize such periodic operators 
by the Bloch-Floquet transform defined next. 

Definition 3 The Bloch-Floquet transform 7 : ^^(Z) ^ L'^^n) is defined by 
{F<PUk) = — ^ '^«+-^ e<'^+'"^)^ n G 0, . . . , iV - 1 , A; G , 



where Tn = {-jj, §] and \Tn\ = ■ Its inverse : L2(TAr) C-^ i^{Z)^C^ is given 
by 

(J^-V)n = /— - / dk (f)nmo<iN{k) e^''''' , neZ. 
a/ITatI 7tjv 

Proposition 11 One has the following properties: 

(i) J-""^ = J-"*, namely T is unitary. 

(ii) Let T 6e the shift on £^(Z) defined by T(j)n = 0n+i and let T^^^ be the cyclic shift on C^. 
Then 

iTTT*(l)Uk) = e-'' (T.,.0)„(A;) , n G Z , A; G . 

Proof. This follows from direct computations. □ 

For any k G Tn and Sj = S{aj, Uj, Vj) G U(2L)inv let us set 



S,{k) - 

Now V^''(/c) and W^^(/c) are defined as Section [8] using Sj{k) instead of Sj and then I[J^''(/c) = 
V^"(/c)W^'(/c) is a finite periodic scattering zipper, the spectrum of which can be calculated by 
the technique of Section |8| By the following result this allows to calculate the spectrum of U*""". 

Theorem 7 The operators V""', and U""" are fibered after Bloch-Floquet transformation, 

precisely 

/•© r® r® 

jryper^r* _ / dkY'j^'{k) , J^W^^J^* = / dkW^'{k) , = / dkWj^'ik). 

Jjm -^TTAr Jtn 

Therefore, 

= U aiWj^^k)) . 
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Proof. By definition, for G ^^(Tjv) C^: 



m£2 



In case n = 2j — 1, one obtains from the structure of V'" 
(^V--J-»2,-i(A:) 



==r ^ (tt2j(-^*0)2j-l+mAr + (1 - a2ia2j)^^2i(^*0)2i+miV 



g«(2j-l+mAr)fc 



= 02J-1 + (1 - «2i«2i)'f^2i e" </>2j 

By the same calculation for n = 2j, 

(^V-'^*0)2,(fc) = V^2,(l-a;,«2,)^e^' 02,-1 - V^2,a;,?72, 02, . 
Together this shows 

jryp-jr* = dk S2{k) ® 54(A;) ® ■ ■ ■ 5^(A;) = dk V^^(A;) . 
The second equality is proved using Proposition [TTj^ii) and previous computation: 

= r dkT^,^{S^{k)®---®SN-i{k))T-l = rdkW-{k). 
As the product of fibered operators is fibered, this also implies the formula for V-""". □ 

Appendix A: Mobius transformations 

This appendix resembles some basic properties of the Mobius transformation as they are used 
in the main text. A lot of references to the literature can be found in |DPS] . The Mobius 
transformation (also called canonical transformation or fractional transformation) is defined by 

T-Z = {AZ + B){CZ + D)-\ r = (^^ G G1(2L,C), Z G Mat(L x L, C) , (43) 

whenever the appearing inverse exists. If T is ^/-unitary and Z G M.^, then T- Z exists and is in 



Ml (see Appendix B for the definition of the upper half-plane Ml). For T as in (43) and as long 
as the appearing inverse exists, the inverse Mobius transformation is defined by 

W.T = {WC - A)-^ {B - WD) , W G Mat(L x L, C) . (44) 

The Mobius transformation is a left action, namely (TT') ■ Z = T ■ {T' ■ Z) as long as all objects 
are well-defined. The inverse Mobius transformation is a right action in the sense of the following 
proposition, the algebraic proof of which is left to the reader. 
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Proposition 12 Under the condition that all the Mohius und inverse Mohius transformations as 
well as matrix inverses below exist, one has the following properties. 

(i) W = T-Z ^ W:T=Z 

(ii) W : {TV) = {W:T):r 

(ill) w ■.r = T-^-w 

Appendix B: Riesz-Herglotz representation theorem 

Let Ml denote the upper half plane of matrices Z G Mat(L, C) such that Qm{Z) = i{Z* — Z) > 0. 
It is well-known {e.g. Section 4.5. of |Sim2| ) that the Cayley transform maps Ml via Mobius 
transformation to the Siegel disc D^, namely C-Ml = D^. An analytic function 2; G Hi i— )• G{z) G 
Ml is called a Herglotz function. Then F{z) = G{C* ■ z) is an analytic function on the unit disc 
Di having positive imaginary part. If, moreover, -F(O) = «1, then such a function is called a 
Caratheodory function. The scalar version of the following classical theorem can be found in text 
books such as [Lax]. The matrix version is an immediate corollary of it. 

Theorem 8 (Riesz-Herglotz representation theorem) Let F : D — )■ Mat(L x L, C) he analytic 
satisfying '^m{F{z)) > and F{0) = il. Then there exists a unique matrix-valued probability 
measure fi on §^ such that 
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